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We study the conductivity and shear viscosity tensors of a strongly coupled M — 4 super- Yang- 
Mills plasma which is kept anisotropic by a S parameter that depends linearly on one of the spatial 
dimensions. Its holographic dual is given by an anisotropic axion-dilaton-gravity background and 
has recently been proposed by Mateos and Trancanelli as a model for the pre-equilibrium stage 
of quark-gluon plasma in heavy-ion collisions. By applying the membrane paradigm which we 
also check by numerical evaluation of Kubo formula and lowest lying quasinormal modes, we find 
that the shear viscosity purely transverse to the direction of anisotropy saturates the holographic 
viscosity bound, whereas longitudinal shear viscosities are smaller, providing the first such example 
not involving higher-derivative theories of gravity and, more importantly, with fully known gauge- 
gravity correspondence. 
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Introduction. Hydrodynamic simulations of heavy- 
ion collisions suggest that the produced quark-gluon 
plasma is behaving like an almost perfect fluid with a 
ratio of shear viscosity over entropy density not far from 
the famous result H/Att associated with the membrane 
paradigm of black holes [2] and which holographic gauge- 
gravity duality maps to the corresponding quantity of 
maximally supersymmetric Yang-Mills theory in the limit 
of infinite color number and infinite 't Hooft coupling 
Oil]. This value has been conjectured to form the lower 
bound for any realistic matter |5j. It was found to be 
saturated universally [SJ [7] in dual theories involving an 
isotropic horizon described by Einstein gravity. Values 
above this bound are obtained when corrections due to 
finite coupling strength are included [S], but it has been 
shown that values violating the bound can arise in higher- 
derivate gravities [9] , although so far no complete gauge- 
gravity correspondence has been established for finite vi- 
olations. 

The possibility of nonuniversal shear viscosity within 
two-derivative gravity has been considered before in grav- 
ity duals of noncommutative gauge theories [lOj and p- 
wave superfiuids fllj, where shear modes in addition to 
the universal purely transverse one exist. For p-wave su- 
perfiuids, nonuniversal values above the bound have been 
recently obtained in [12]. 

Anisotropics are of particular importance to the 
nonequilibrium dynamics of the quark gluon plasma in 
heavy-ion collisions. At weak coupling, they give rise to 
non-Abelian plasma instabilities |13j and turbulent be- 
havior [2] that has been argued to effectively reduce 
the shear viscosity [T5], which is parametrically large in 
a weakly coupled plasma [T^]. More directly accessible 
experimental signatures such as anisotropic photon and 
dilepton emission [17], momentum broadening of jets [18] 
as well as quarkonium dissociation [19j have been studied 



in the hard-loop effective theory "20] for a weakly coupled 
anisotropic plasma where on sufficiently small time scales 
the anisotropy can be taken as stationary. 

While holographic models of the evolution of strongly 
coupled quark-gluon plasma towards isotropy are of im- 
mediate interest and have been studied extensively [21j . 
further insight may be obtained by models where the 
anisotropy is temporarily fixed. The first such attempt 
was presented by Janik and Witaszczyk in [52] , where the 
gravity dual involves a naked singularity. In |23j we have 
studied electromagnetic signatures of this model which 
are qualitatively similar to weak-coupling results at high 
frequencies, but which involve vanishing conductivities in 
the limit of zero frequency and absence of hydrodynamic 
behavior. 

Anisotropic axion-dilaton-gravity solution. Recently, 
Mateos and Trancanelli [24l [25] have adapted the spa- 
tially anisotropic duals of Lifshitz-like fixed points of [55] 
to a setup with AdS boundary conditions and finite tem- 
perature. This provides an anisotropic version of an 
TV = 4 super- Yang-Mills plasma where the anisotropy 
is kept stationary by deforming the gauge theory by 

SS=^ J e{z)Tr FAF (!) 

with 6{z) — 2'Kaz depending linearly on one of the spa- 
tial dimensions. The constant a is related to the density 
of D7-branes (homogeneously distributed along z) that 
are dissolved in the bulk of the dual theory (they do not 
extend to the holographic boundary) . The dual bulk the- 
ory is effectively 5-dimensional axion-dilaton gravity with 
negative cosmological constant and bulk action 

^^"'^=2;^^/^^^ + ^'- — 2 J 

where — 8ttG — Att'^/N^ and axion x — clz. Here we 
have dropped a factor of the ten-dimensional solution 



2 



Ai X and set its radius L = 1. The remaining five- 
dimensional geometry which asymptotes to AdSs for u — )■ 
can be given by a Une element of the form 



-0(«)/2 



- F{u)B{u) dt 
+ dy^ +niu)dz^ 



du" 



(3) 



The functions B, T and Ti. can be written in terms of the 
dilaton (/>, which itself has to satisfy a third-order non- 
linear differential equation in u pS) . For nonvanishing 
anisotropy parameter a there is an anisotropic horizon 
J-{uh) — with H > 1, which is completely regular. This 
will allow us to apply the membrane paradigm along the 
lines of 7 for studying conductivities and shear viscosi- 
ties, both of which have two independent components at 
nonzero anisotropy. 

The thermodynamics of this setup and its instabilities 
are discussed in |25j . Upon holographic renormalization 
which brings in a reference scale fj,, the stress tensor of 
the gauge theory has the form 



=diag(e,Pi,Pi,F,), 



(4) 



with a conformal anomaly (T^) oc a"" 
sity is given by 

s = {e + P^)/T 



The entropy den- 



(5) 



with T = IJ'' {uh)\^/Bh / 4:Tr . Unlike the strictly positive s 
which continues to be given by a quarter of the horizon 
area over spatial volume, e and the pressure components 
depend separately on T/fj, and a//x and can become neg- 
ative in various regions of the phase diagram. When \a\ 
is increased from the isotropic limit a — 0, the pressure 
anisotropy first becomes oblate (P^ < Pj_), while prolate 
anisotropies can be obtained for larger |a|. For small a 
there is a qualitative similarity to the plasma instabilities 
in an anisotropic weakly coupled plasma in that the ho- 
mogeneous phase is unstable against filamentation of the 
distribution of the dissolved D7 branes along z [IS] . For 
sufficiently large a the instabilities against filamentation 
however disappear, and also various components of the 
stress tensor (but not the entropy) can become negative. 

Conductivities. Including a standard Maxwell action 
for a U(l) gauge field in the bulk with coupling constant 
(7eff leads to a conjugate momentum with respect to evo- 
lution in the coordinate u that defines the current 



dC 



(6) 



As in [7], regularity at the horizon in terms of the 
Eddington-Finkelstein coordinate v for infalling ob- 
servers defined by 



dv = dt — 



-du 



-gtt 



(7) 
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FIG. 1: DC conductivities along and transverse to the direc- 
tion of anisotropy as a function of the anisotropy parameter 
a/T. 



relates linearly to the electric field Fu, so that the 
horizon acts as a conductor. The anisotropic metric ([s]) 
however leads to two different conductivities 



= .9off 



(8) 
(9) 



where 7 = g/guduu- In the limit of vanishing frequency 
uj and wave vector q we have d^j^ — >■ which allows us 
to conclude that az,± give the DC conductivities of the 
boundary theory. The numerical results as a function of 
the anisotropy parameter a/T are given in Fig. [I] 

Because az^± are determined by geometric data, they 
do not depend on fi and thus only on the ratio a/T. 
Whether a given value a/T corresponds to oblate or pro- 
late anisotropy however depends on T/ fj, |25] . Our finding 
that az/cr^ is always smaller than one may seem surpris- 
ing, but we have in fact found an extreme version of this 
phenomenon in the singular geometry of |22j . While in 
this case both DC conductivities vanish, the ratio of AC 
conductivities CTz (uj) / a± (uj) — >■ as w — J' for both oblate 
and prolate anisotropies (see Fig. 11 of [15]). 

Following the procedure of [27. , one can also obtain two 
related diffusion constants Dz^± = <Jz,±'^~^ with com- 
mon charge susceptibility 



du- 
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(10) 



We have checked [55] the results for Dz.i_ by numerical 
calculation of the lowest quasi-normal mode in the diffu- 
sive channel, thus verifying that the membrane paradigm 
is indeed valid for conductivities of an anisotropic hori- 
zon. 

Shear viscosities. In an anisotropic fiuid with axial 
symmetry, the viscosity tensor |29j rjiji^i has five inde- 
pendent components, two of which are shear viscosities 
which we denote as 



(11) 
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and which can be extracted by Kubo formulas from the 
corresponding retarded correlation functions of the stress 
tensor. In the gravity dual we need to consider metric 
fluctuations /i^^, for which we adopt the gauge h^u — 
0. Expanding the action ^ to second order in /i^^ we 
obtain 



1 



167rG 



-(0) 
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where we introduced 



yxx 



(12) 



(13) 



Considering first Fourier modes g) = hy{u,q) 

and their conjugate momenta Il±(u,q) we find that re- 
quiring regularity at the horizon as above for vector per- 
turbations gives the purely transverse shear viscosity 



V-L = 



iujTl>j_{uh,q) 



167rG 



s 

An 



(14) 



which reproduces the universal value for Einstein grav- 
ity with isotropic horizons. Because the radial evolution 
becomes trivial as (7'^ — > 0, 

d^IlAu,q) = ^{g*'uJ^+g^-^ql)^/j^{u,0) ^ 0, (15) 

the (momentum independent) result of the membrane 
paradigm gives the purely transverse shear viscosity of 
the boundary theory at = 0, exactly as in the univer- 
sality proof of [7j. 

Turning to nontransverse metric fluctuations we find 
that the radial flow is still trivial in the limit — > for 
the case of = 



"hzx because 



duTlL{u,q) = 



-9 9x 



IGttG 



'-(g''u:^+gyyql)i,L{u,G)-^0 



for any a (while 



(16) 

would have a nontrivial flow oc 



a ). We can therefore use the membrane paradigm to 
calculate the longitudinal shear viscosity from 

^\\{uh,q) \/l{uh) 9xx{uh) 



iuj^j\\{uh,q) 
s s 



167rG gzz(uh) 



(17) 



for a ^ 0, which is plotted in Fig.[2]as a function of a/T. 
Note that the ratio ri\\/rij_ turns out to be identical to 



We have also checked (17) by a direct numerical eval- 
uation of the Kubo formula 



lim lim Im — — '—^ . 



(18) 
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FIG. 2: Transverse and longitudinal shear viscosities over 
s/47r as a function of the anisotropy parameter a/T. 



Momentum diffusion. As a further corroboration, we 
have also studied the diffusive behavior of the relevant 
shear modes. The longitudinal shear viscosity is re- 
lated to the diffusive pole of the shear mode Z3 = 
g^^ (^qhtx{t, z,u) + ujhzx{t, z,u)'j at small but nonvanish- 
ing momenta, where q is oriented in the z-direction in 
order to preserve 0(2) symmetry which avoids a coupling 
to axion modes (the general case will be discussed else- 
where [55] ) . We can use the trick of [23 of Kaluza-Klein 
dimensional reduction in the a;-direction, which maps the 
problem to that of U(l) gauge fields with a u-dependent 
effective coupling. This yields the diffusion constant 



(19) 



Ufi) Jo 



du 



i Bju) 

n{u)' 



We have also checked this result by calculating numeri- 
cally the lowest lying quasinormal mode of the retarded 
propagator of the shear mode Z^iu) which satisfies the 
ordinary differential equation 



{9tt9'zz - 9zz9'tt) 



(q^9tt 



^^9z 



9uu 



9tt9zz {q'^gtt+uj'^9zz) 



+ 



Ag 



9x 







^3 ~ 9uu 


\9tt 


9zzJ 



z. = 0. 



(20) 



We can now compare ( 19 ) with the analogous cal- 
culation of the shear mode Zi = g^^ (^qhty{t,x,u) + 
ujhty{t,x,u)^ with momentum oriented in the x direc- 
tion, yielding (the breaking of SO (2) invariance in this 
case does not lead to complications from axion modes) 



D 



^v/HK) 



xy 



, 3 3*W B{u) 
du u-^e 4 \ / 

V 



VI- 



47rT 



(21) 



4 



where the last equahty follows from ( 14 ) together with 
([5]), but is also reproduced numerically. Taking the ratio 
of Eqs. (19) and (21) leads again to the result 



< 



i-n-niuh) 47r 



for a 7^ 0. 



(22) 



Conclusion. As mentioned above, a result for r/|| that 
deviates from the universal result of Einstein gravity with 
isotropic horizons has been recently obtained in a dual 
model for p-wave superfluids [T^], however one that is 



above s/in. To our knowledge, the result (22 ) is the first 
example of a shear viscosity that falls below this value 
without recourse to higher-derivative gravity and more- 
over with fully known gauge-gravity correspondence. 

Parts of the phase diagram in a/jj, and T/ ^ of this the- 
ory exhibit instabilities against filamentation into inho- 
mogeneous mixed configurations with aiso = and some 
flani > a, while for sufficiently large a the anisotropic 
homogeneous phase is thermodynamically stable |25j . 
Strictly speaking, the above calculation of hydrodynami- 
cal transport coefficients is limited to stable homogeneous 
phases, but we have numerically verified their existence 
on lines of constant a/T ^ 1.3. 

It is not clear what implications, if any, our result of 
^11 /'/-L < 1 ioT a strongly coupled stationary anisotropic 
plasma would have for real nonequilibrium evolution. 
Intrinsically anisotropic hydrodynamic descriptions of 
quark-gluon plasma have recently been proposed in [50] . 
however without allowing for different shear viscosity 
components. It would be interesting to study the ef- 
fects of different shear viscosities in a hydrodynamical 
evolution, and also what ratios Tyy/Tyj. a weakly coupled 
anisotropic non-Abelian plasma would give rise to. 
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